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Abstract 

This paper studies the problem of pole assignment for symmetric and Hamiltonian 
transfer functions. A necessary and sufficient condition for pole assignment by complex 
symmetric output feedback transformations is given. Moreover, in the case where the 
McMillan degree coincides with the number of parameters appearing in the symmetric 
feedback transformations, we derive an explicit combinatorial formula for the number 
of pole assigning symmetric feedback gains. The proof uses intersection theory in pro- 
jective space as well as a formula for the degree of the complex Lagrangian Grassmann 
manifold. 
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1 Introduction 



One of the best known inverse eigenvalue problems from linear system theory is that of pole 
assignment, i.e. to find a static output feedback gain for a given linear system such that 
the closed loop poles of the system coincide with a specified subset of the complex plane. 
Moreover, in the case of finitely many solutions, a formula for the number of pole assign- 
ing feedback transformations is desirable. Early contributions on the subject were obtained 
by e.g. Davison and Wang [2j and Kimura [20], who derived sufficient conditions for the 
solvability. However these conditons were far from being necessary as well. In a series of 
pioneering papers fHl 121]; R- Hermann and C. F. Martin applied tools from algebraic 
geometry to obtain necessary and sufficient conditions, valid for a generic class of systems 
and for complex feedback transformations. Their approach was based on the dominant mor- 
phism theorem [Chapter AG, §17, Theorem 17.3] |2j from complex algebraic geometry. A 
second breakthrough was subsequently made by R. W. Brockett and C. I. Byrnes [S], who 
used intersection theoretic arguments and the Schubert calculus on Grassmann manifolds 
to count the number of pole assigning complex feedback transformations. By refining these 
algebraic-geometric approaches of Hermann and Martin, and Brockett and Byrnes, a num- 
ber of fundamental contributions on the subject were made that finally led to a solution of 
the problem in the real case, with important contributions due to jHl El I2H1 ESj. For an 
excellent survey paper on this subject, written from a control-theoretic point of view, see 
e.g. |3]. More recently, various intersection theoretic tasks related to the Schubert calculus 
have been studied in the algebraic geometry literature; see e.g. [121 dH IS2]- The focus of 
most of the investigations has been so far on the unstructured case, where no underlying 
symmetries for the involved transfer function or for the associated feedback transformations 
are imposed. However, transfer functions with symmetries occur naturally in various appli- 
cation areas, such as in network theory or mechanics. For example, the transfer functions 
G{s) of linear RLC - circuits, consisting solely of restistors, capacitors and inductive ele- 
ments are symmetric, i.e. they satisfy G(s)* = G{s). In mechanics, the transfer functions 
of linear Hamiltonian systems are characterized by the symmetry relation G(— s)* = G{s), 
while second order mechanical systems of the form 

Mx = Nx + Bu, y = B^x 

yield symmetric Hamiltonian transfer functions, satisfying 

G{s)=H{s'),H{s) = H{sy; 

see e.g. [H [HI IHl IH]- For such structured systems it is reasonable to restrict the class of 
admissible feedback transformations to those that preserve the symmetry properties of the 
transfer functions. Therefore the known results on pole placement on unstructured systems 
do not apply in these cases and require instead a new approach. 

In this paper we start an investigation of the pole placement problem for n x n symmet- 
ric transfer functions G{s) = G'(s)*, arising in electrical network theory, and Hamiltonian 
transfer functions. For both types of systems the natural class of admissible output feedback 
tranformations are the symmetric ones F = F^, yielding a symmetric closed loop transfer 
function 

Gpis) := {In-Gis)F)-'Gis). 
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As the number of free parameters occuring in the symmetric feedback matrices F is n{n + 
l)/2, a necessary condition for generic solvabihty of this output feedback problem is that 
the McMillan degree 6 of the transfer function G satisfies 6 > ("'^^) in the symmetric case, 
and 6 > n{n + 1) in the Hamiltonian case. In fact, we show that generically for complex 
symmetric output feedback transformations this condition is also sufficient. Moreover, for 
the limit case 6 = ("2^) (or 6 = n(n + l)), we derive an explicit combinatorial formula for the 
number of complex symmetric output feedback gains that place the poles at given points. 
Our formula coincides with that of the degree for the complex Lagrangian manifold, given 
in [Si. 

In the real case such complete results can not be expected. In fact, the symmetry of the 
transfer functions then imposes a priori limitations for the possible pole locations of such 
systems. This has been observed in [22]; where it is shown for symmetric transfer functions 
that - in the special case that the Cauchy index of G coincides with the McMillan degree - 
then generically real symmetric output feedback pole assignability holds if and only n > S. 
Of course, in most applications we have n < S and therefore the description of the set of 
poles that can be achieved by real symmetric output feedback becomes a complicated and 
nontrivial task. 

2 Complex symmetric and Hamiltonian realizations 

In this section we recall some basic facts concerning complex symmetric and Hamiltonian 
transfer functions, respectively and associated signature symmetric and Hamiltonian real- 
izations. Let C denote the field of complex numbers. A complex rational transfer function 
G{s) G C(s)"^" of McMillan degree 6 is called symmetric, or Hamiltonian, respectively, if 

G{s) = G{sY, or G{s) = G{-s)\ respectively, 

holds for all s G C. A complex symmetric realization is a linear system of the form 

X = Ax + Bu, y = B^x, 

where A G C^^"^ is symmetric, i.e. A* = A, and B G C^^". Similarly, a Hamiltonian 
realization is a linear system 

X = Ax + Bu, y = Gx, 
where A G C^^^ B G C^^", G G C"^'' satisfies 

AJ = [Ajy, C* = JB 

and 

I 
-I 

denotes the standard symplectic form on C^^"^. In particular, Hamiltonian systems have 
always even McMillan degree 5. 

Complex symmetric realizations are the natural class of realizations for complex sym- 
metric transfer functions. In fact, they are the proper analogue of signature symmetric 
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realizations of real rational transfer functions, appearing in network theory. Over M, real 
symmetric realizations correspond to linear models of RC— networks, constructed entirely 
using capacitors and resistors. The real symmetric transfer functions defined by them are 
characterized by the property that the Cauchy-Maslov index coincides with the Mcmillan 
degree, [Hin]- 

The following variant of the Kalman realization theorem is well-known; see e.g. pi llUj. 
Recall that the complex orthogonal group 0(5, C) is the matrix group consisting of all com- 
plex 5x5 matrices satisfying 5*5** = Is- Given any complex realization (A, 5, C) of a 
symmetric transfer function G{s) = C{sl — A)^^B, note that (A*, C*, S*) is also a realization. 

Proposition 2.1. Let G{s) = G{sy be an n x n strictly proper, complex rational transfer 
function of McMillan degree 6. Then 

(1) G{s) has a controllable and observable complex symmetric realization (A, B, C) = 
{A\G\B'). 

(2) If {Ai, Bi,Ci) , i=l,2, are two controllable and observable complex symmetric realiza- 
tions of G{s), then there exists a unique complex orthogonal transformations S G 
0(5, C) such that (^2, 52,02) = {SAiS-\SBi,CiS-^). 

In the literature usually only the real case of the above result is proven, where the state- 
ment is actually slightly different due to the presence of signature symmetric realizations. 
In the complex case the result simplifies to the one given here. For the sake of completeness 
we include the proof; see also [TU]. 

Proof. If {A, B, O) is a minimal realization of G{s) then, by symmetry of G, also (A*, O*, 5*) 
is a minimal realization. Applying Kalman's realization theorem implies the existence of a 
unique invertible complex 5x5 matrice S with 

{A\G\B') = {SAS-\SB,GS-^). 

By transposing this equation and using the uniqueness of S we conclude that S = S^. It 
is a well known fact from linear algebra that every complex symmetric invertible matrix 
has a representation 5* = XX* by a complex invertible matrix X. Moreover, X is uniquely 
determined up to right factors XT where T G 0(5, C). Then {XAX~^,XB,GX~^) is a 
complex symmetric realization, which completes the proof. □ 

There is a similar realization theorem for Hamiltonian systems, for which we refer to the 
literature; see e.g. [HI El- Static linear output feedback can be meaningfully defined for such 
systems only through symmetric gain matrices. Thus an output feedback transformation 

u = Fy + V 

with the closed system 

x = {A + BFB^)x + Bu,y = B^x 

preserves the complex symmetry of the realizations if and only if F = F*. Thus we define 
two complex symmetric realizations {Ai, Bi, Gi) to be symmetric output feedback equivalent 
if and only if there exist 5 G 0(5, C), F = F* G C"^" with 

(A2, B2, O2) = {S{Ar + BrFB\)S-\ SB^, G^S-'). 
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Equivalently, if and only if for the associated transfer functions Gi{s): 

G2{s) ■.= iIn-G^is)F)-'G,is). 
Similarly, output feedback for Hamiltonian systems 

x = {A + BFC)x + Bu,y = Cx 

preserves the Hamiltonian properties of the realization if and only if F = F*. Thus in both 
cases we have to focus on symmetric output feedback. 

We note some elementary geometric properties of the set of complex symmetric transfer 
functions that will be important in the subsequent development; see e.g. jHj for providing 
some of the details for the proof of the subsequent theorem. We omit a full proof as it would 
take us to far apart from the subject. 

Proposition 2.2. Let SRats,n{'C) and Hams^ni'C) , respectively denote the sets of strictly 
proper, complex symmetric and Hamiltonian, respectively, nxn transfer functions of McMil- 
lan degree 6. Then SRats^ni'^) one/ Ham^ni^), respectively, is a smooth complex manifold of 
complex dimension 5{n+l), and dimension 6n respectively. Moreover, they are nonsingular 
irreducible quasi-affine varieties. 

In particular, there is a canonical notion of "genericity" for complex symmetric or Hamil- 
tonian transfer functions. Explicitely, a property E of complex symmetric transfer functions 
is called generic, if the set defined by E 

{G E SRats,n{'C) I G has property E} 

is a Zariski-open subset of SRats^miC). Equivalently, this can be also expressed in terms of 
complex symmetric realizations. 

3 Main result 

After these preliminaries we can now rigorously formulate and proof the main technical 
results of this paper. Let G{s) be an n x n complex symmetric or Hamiltonian transfer 
function, i.e. = G{s) or G(— s)* = G{s), respectively. Assume that G{s) is strictly 

proper and has McMillan degree 6. The complex symmetric eigenvalue assignment problem 
then asks the following question: 

Problem 3.1. Given an arbitrary monic polynomial if{s) G C[s] of degree 6 (f{s) = <f{—s) 
is assumed to be even in the Hamiltonian case). Is there an nxn complex symmetric matrix 
F such that the closed loop transfer function 

GF{s):={In-G{s)Fr'G{s) 

has characteristic polynomial ip{s) , i.e. the poles of Gf{s) are the zeroes ofip{s)? 
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If for a particular symmetric (Hamiltonian) transfer function G{s) Problem 13.11 has an 
affirmative answer we will say that G{s) is pole assignable in the class of complex symmetric 
(Hamiltonian) feedback compensators. We say that G{s) is generically pole assignable, if the 
problem is solvable for a generic choice of admissible polynomials f{s). 

Similar to the situation of the static pole placement problem |3l ES] and the dynamic 
pole placement problem j2Z], Problem 13. II turns out to be highly nonlinear and techniques 
from algebraic geometry will be required to study the problem. The first main result is in 
the spirit of Hermann and Martin, by deriving a generic necessary and sufficient condition 
via the dominant morphism theorem. 

We prove some lemmas first. Let tt{A) = (an, . . . , ass) be the projection onto the diagonal 
entries of an 5 x 5 matrix A. In the sequel we will identify C'' with the complex vector space 
of row vectors. For any symmetric matrix L, define 9^ : 0{6, C) through 

e^{S) = TciSLS^^). 

As 0(5, C) is a Lie group, its tangent space at the identity matrix / is given by the Lie 
algebra of complex skew-symmetric matrices 

so{6, C) = {X e C^^^ I X + X* = 0}. 

Moreover, the Jacobian d6f of 6^ at / is given by 

def : so{5, C) ^ V, def{X) = tt{XL - LX), 

where 

5 

y = {(xi,...,X5)GC^l5^x, = 0}. 

1 

For any 5x5 matrix L, the graph Q{L) of L is defined as a graph with 5 vertices such 
that there is a path from vertex i to vertex j if and only if the ijih entry of L is none zero. 
It is a well-known fact from linear algebra, that the graph Q{L) is connected if and only if 
L is irreducible, i.e. if and only if there exists no permutation matrix P such that PLP"^ 
is block diagonal. We use this fact together with an idea developed in Lemma 2.5] to 
prove the following equivalent characterization. 

Lemma 3.2. The Jacobian d9f is surjective if and only if the associated graph G{L) is 
connected. 

Proof. By inspection, the derivative d6f is not surjective if and only if there exists a nonzero 
diagonal matrix Z of trace zero, such that for all X G so{6, C) 

trace(Z(XL - LX)) = trace((LZ - ZL)X) = 0. 

By symmetry of L,Z we have LZ — ZL G so{6,C). Since the trace function defines a 
nondegenerate bilinear form on so(5, C), the condition iia,ce{{LZ — ZL)X) = is equivalent 
to LZ = ZL. Since Z is a nonzero diagonal matrix of trace zero, there is a permutation 
matrix P such that Z := PZP^^ = block diag(ai/i, . . . , a^/fc) with k > 2 and a^'s distinct. 
Let L = PLP~^. Then LZ = ZL is equivalent to LZ = ZL, which is equivalent to L 
being block diagonal. But from the above remark this is equivalent to the graph Q (L) being 
disconnected. The result follows. □ 
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Lemma 3.3. Let L be a nonzero complex symmetric matrix such that 7r(L) = 0. Then 
there is a family of orthogonal matrices S{e) e 0{S, C), e > 0, with S{0) = I such that the 

matrix L{e) := S{e)LS{e)^^ has the properties that 7r(L(e)) = and dip^^"^^ is surjective for 
alle e (0,7r/2). 

Proof. If G{L) is connected, then by the previous lemma the choice S{e) :— I does the job. 
Thus it suffices to prove that Q{L) not connected imphes that then we can find a family of 
transformations >S'(e), such that TT{S{e)LS{e)^^) = and the largest connected subgraph of 
Q{S{e)LS{e)^^) contains more vertices than that of Q{L) for all < e < n/2. 

Note that 7r(L) =0 and L imply that the largest connected subgraph of ^(L) must 
contain at least 2 vertices. Assume that the largest connected subgraph of G{L) contains k 
vertices, 2 < k < S. Without loss of generahty, assume 



Li 
L2 



where the graph of the k x k sub-matrix Li is connected. Write 



Ln a 
a* 



and L2 



/3* 
P L22 



where Lu and L22 are sizes {k — 1) x {k — 1) and {6 — k — 1) x {S — k — 1) respectively. By 
irreducibility of Li we have a 7^ 0. Thus L has the form 



a* 





a 





p 



Let 



-^(fe-l)x(fe-l) 









cose 
sin e 







L22 


— sin e 
cos e 










Then 



S{€)LS-\€) 



[cos e)a 
(sin e)a^ 




(sine) a 

(-sine)/?* 
(cose)/3* 
L22 



[cos e)a 



[—sme)P (cose)/3 

For the graph Q{S{e)LS~^{e)) with < e < tt/2, the vertices {1, . . . ,k} are still connected 
and the vertex A; + 1 is symmetrically connected to at least one of the first k vertices. Thus 
the vertices {1, . . . , /c + 1} are connected. □ 

Lemma 3.4. Let £ be a linear subspace of complex symmetric matrices of dimension 5, and 
C ^ s/((5, C). Then there exists an orthogonal matrix S e 0(5, C) such that tt \scs-i is one 
to one, and onto. 
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Proof. The proof goes by recursively constructing a basis {Li, . . . , Ls} of C such that 
{ni^SLiS^^), . . . ,tt{SLsS~^)} are hnearly independent for a suitable complex orthogonal 
matrix S G 0(5, C). First note, that we can modify any basis of C into a basis C^^^ : = 
{Li, . . . , Ls} of £ such that Li ^ sl{6, C), and Li e C), for i = 2, . . . , 5. In fact, if 
{Ki,...,Ks} denotes any basis of £ with trace(i^'i) ^ 0, then {Li := Ki,L2 := — 
C2K1, . . . ,Ls := -ft'^ — c^i^'i}, Cj := trace(-ft'j)/ trace(_ft'i), is as desired. By construction of 
Li, then 7r(Li) 7^ 0. 

Let {Li, . . . , L^} be a basis of £ such that Li ^ C), and Lj e s/(5, C), for i = 2, . . . , 5. 
Then dim span{7r(Li), . . . , irlLs)} := k > 1. U k < S, then by re-ordering the indices we can 
assume that {7r(Li), . . . , 7i{Lk)} are linearly independent, and 

k 

n{Lj) = ^ Cij7i{Li) for j = /c + 1, . . . , 5. 

1=1 

By replacing Lj with Lj — Yli=i^ij^i '^^^ further assume that Tx{Lj) = for j = k + 
1,...,6. It is thus sufficient to show that if there is an orthogonal matrix S such that 
the matrices {Mj := SLjS~^,j = 1, . . . ,6} have the property that {7r(Mi), . . . , 7r(Mfc)} are 
linearly independent, and 7r(Mj) = 0, j = k + 1, . . . ,6, for some k < n, then we can find an 
orthogonal S such that 

{7c{SM,S-'), TiiSM.S-'), 7riSM,+,S-')} 

are linearly independent. 

By Lemma 13.31 there exists G 0{S,C) arbitrarily close to the identity matrix such 
that TT{SeMiS~^), . . . ,7!'{SeMkS~^) are linearly independent and the graph Q{SeMk+iS~^) is 
connected. By replacing Mj with SeMiS~^ we can assume further that dO^^^^^ is onto V . 
Then there exists a skew-symmetric matrix X such that 

7r(XMfc+i - Mfe+iX) ^ span{7r(Mi), . . . , 7r(Mfc)}. 

Let 

S{e) = exp(eX). 

Then S{e) is orthogonal for all e, and 

S{e) = I + tX + higher order terms. 
The Taylor series expansions of {T[{S{t)MiS{t)^^)} have the forms 

Ti{S{t)MiS{t)-') = 7r(M,) + A(e), ^ = 1, . . . , fc, 

and 

7r(5(e)A4+i5(e)-i) = e (7r(XA4+i - M^+iX) + /3fe+i(e)) 

where /3j(e) are continuous with respect to e and /5j(e) ^ as e ^ for i = 1, . . . , A; + L 
Since {7r(Mi), . . . , 7r(Mfc), t^^XM^^i — M^+iX)} are linearly independent, for sufficient small 
e > 0, {7r(Mi) + A(e), . . . , 7r(Mfc) + /5fc(e), 7r(XMfe+i - M^+iX) + /9fc+i(e)} are also linearly 
independent, i.e. {7r(S'(e)MiS'(e)~"^, . . . , 7r(S'(e)MfcS'(e)"\ 7r(S'(e)Mfc+iS'(e)^^} are linearly in- 
dependent. □ 
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Theorem 3.5. IfG{s) is a symmetric (or Hamiltonian) transfer function of McMillan degree 
5 > (^2^) (or 6 > n(n + 1) ), then G{s) is not pole assignable in the class of (real or) complex 
symmetric feedback compensators. 

When 6 < ("2^) (or 6 < n{n + 1)), then there is a generic set of n x n symmetric (or 
Hamiltonian) transfer functions of degree S which are generically pole assignable via complex 
symmetric feedback compensators. 

Proof. We only give a sketch of the proof, as the arguments based on the dominant morphism 
theorem are well known from pTHl2Hj. Note, however, that there is serious gap in the proof of 
[SB] for the pole placement result on Hamiltonian systems because it is not proved that the 
set of generically pole assignable Hamiltonian systems is non empty. In fact, a construction 
of such an example is not completely trivial and depends on our previous lemmas. 

The first claim follows immediately from a standard dimension argument, as the vector 
space Sym(n) of complex nxn symmetric matrices has dimension ("^ ) . For the second claim 
we note that the set of generically pole assignable systems is a Zariski open subset of the 
nonsingular, irreducible quasi-affine variety of symmetric or Hamiltonian transfer functions, 
respectively. Therefore we only need to show that this Zariski open subset is nonempty. By 
the Dominant Morphism Theorem, it suffices to find one system whose Jacobian of the pole 
placement map at one point is onto. 

Note, by the Newton formula, that the coefficients of the characteristic polynomial 
det(s/ — A) = + as-is^~^ + ■ ■ ■ ais + ao are related to the traces of powers of A as 
follows: 



as-i = — trace(74) 

as-2 = — -( trace(A^) + trace(y4)) 



ao = —-{tiace^A^) + as-itYa.ce{A^ ^) + ■ ■ ■ + ai trace(A)). 


Therefore for the case of symmetric transfer functions, the pole placement map is equiv- 
alent to the map 



(3.1: 



(p : Sym(?7,) — > 

F I — > ( trace(A + BFB*), ■ ■ ■ , trace(A + BFB^) 

and its Jacobian at is given by 

rf0o(F) = (trace(5F5*),2trace(AfiFfi*),...,5trace(A'^"^fiF5*). 

For the case of Hamiltonian transfer functions, since JAJ = A* and = —I, we have 
{—1)''^^JA^J = {A'^y for /c = 1, 2, . . . , which implies that the characteristic polynomial of 
A is even and 

trace(A'') = holds for all odd A;'s. 
Therefore the pole placement map is equivalent to the map 

^: SymH 

F I — > ( trace(A + BFCf, trace(A + BFC)^, ■ ■ ■ , trace(A + BFCf) ^ ' 
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and its Jacobian at is given by 



dipo{F) = (2trace(ASFC),4trace(A'^5FC) . . . ,6tTa.ce{A^-^BFC). 

We first consider the case of symmetric transfer functions. Let B be any real nonzero 
matrix and C = {BFB^ \ F G Sym(ri)}. Then C (f_ s/((5, C) and dim£ > b. By Lemma IX 
there exists an orthogonal matrix S G 0(5, C) such that ttI^^^-i is surjective. Let D 
diag(l, 2, . . . , 5) and A = S~^DS. Then 



d(\)Q{F) = (trace(5F5*),2trace(A5F5*),...,(5trace(A^-iSFS*)) 
= ( tracelsBFB'S-^), 2 tTace{DSBFB'S 
= Ti{SBFB'S-^)V 



t Q-l\ 



6tTeice{D^-^SBFB'S-^)) 



where 



V 



1 1 
1 2 

1 6 



1 

2S-1 



6 



5-1 



1 
2 





Since ti\scs-'^ is surjective and V is nonsingular, d(j)o is onto. 
For the case of Hamiltonian transfer functions, Let 








B 







and C = [ ] 



where Bi is any real nonzero | x n matrix, and C = {BiFB\ \ F E Sym(?2)}. Then 
C (t sl{6/2,C). By Lemma EiH there exists an orthogonal matrix Si G 0(5/2, C) such that 
TV : SiCS^^ 1-^ Ca is surjective. Let Di = diag(l, 2, . . . , 5/2), 



S 



Si 
Si 



, and D 



Di 
Di 



and A = S ^DS. Note that D, S are Hamiltonian and symplectzic matrices, respectively. 
In particular, A is Hamiltonian. Then 



dMF) 



{2tmce{ABFC),4:tTace{A^BFC), . . . ,StTace{A^-^BFC)) 
{2tracelDSBFCS-^),4:tTace{D^SBFCS-^), . . . ,5tTaceiD^-^SBFCS-^)) 



{2tr&ce{DiSiBiFB\S^^),AtTiice{DlSiBiFB{S 
n{SiBiFBiSi^)U 



t C-1-] 



1 7' 



..,6tTace{Dl-^SiBiFB{S 



where 



U 



1 1 

2 2^ 



s\3 



(f) 



1 

25-1 



(f) 



5-1 



2 
4 










Since ^rl^^^^^^t^-i is surjective and U is nonsingular, dipo is onto. 



□ 
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The second main theorem in this paper deals with the limit case 5 = ("^^), where we 
can prove a more precise statement. 

Theorem 3.6. Let 5 = {^2^) in the symmetric case, and 6 = n(n + 1) for Hamiltonian 
systems. Then for a generic set of n x n symmetric (or Hamiltonian) transfer functions of 
degree 6 the number of pole assigning complex symmetric feedback compensators is finite and 
when counted with multiplicities there are exactly 



din) := 2UJ ^ ' ' = . (3.3) 

1! 3! ■■■ {2n-l)\ nr=o (2^ + 1)"" 

many symmetric compensators as solution. 

One immediately computes d{l) = 1, d{2) = 2, d(3) = 2^ rf(4) = 3 ■ 2^ d{5) = 11 ■ 13 ■ 2" 
and d{6) = 13 ■ 17 ■ 19 ■ 2^*^. The integer sequence d{n) is sequence A005118 in Sloane's 
data bank of integer sequences [3T]. The sequence has several combinatorial and geometric 
interpretations. For the context of this paper it will be important that d{n) is equal to the 
the degree of the Lagrangian Grassmannian, the projective variety of all maximal isotropic 
subspaces in a complex vector space of dimension 2n and this has been recently established 
by Totaro 

As it can be seen from this sequence, d{n) appears always to be even, except for n = 1. 
This is related to the fact, that the symmetric output feedback pole placement problem is 
not generically solvable over the reals. Actually more is true. The sequence 

d{n) := d{n)2'(^) 

is the degree of the spinor variety, the complex projective variety S0{2n + 1)/U{n) J7]; in 
particular d{n) represents an integer sequence again. The sequence d{n) appears under the 
number A003121 in Sloane's data bank |31j . 

The proof of Theorem 13 . 61 will occupy the rest of this section. The proof will necessitate 
a geometric reformulation and several technical lemmas. 

First we will describe the closed loop characteristic equation in a slightly more convenient 
way. Consider a left coprime factorization D~'^{s)N{s) = G{s) of the symmetric or Hamilto- 
nian transfer function G{s). Let F G Sym(n) be an n x n complex symmetric matrix. When 
the feedback law y = —Fu + f is applied then up to a constant factor the characteristic 
polynomial ip{s) is also equal to 



det 



D{s) N{s) 
F L 



(3.4) 



The vector space Sym(n) describing the set of n x n complex symmetric matrices is 
not very well suited to invoke strong theorems from algebraic geometry and intersection 
theory [TT], as these usually require compactness assumptions on the underlying spaces. A 
similar difficulty exists for the static output pole placement problem. Brockett and Byrnes 
showed in |3] how to translate the static pole placement problem into a geometric prob- 
lem. This then resulted into an intersection problem on a compact Grassmann variety and 
methods from classical Schubert calculus 123 ESI could be invoked. 
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We will follow this compactification strategy for Problem 13.11 as well. This will lead 
us to an intersection problem on some projective variety. In order to do so we therefore 
need a good compactification of Sym(n). For this identify the rowspan rowsp [F In] of any 
symmetric matrix F with an element of the Grassmann variety Grass(r;,, C^"). Using the 
Pliicker embedding 

Grass(n, C^") — > P (a^C^'^) = P^, = (^^^ - 1 

we can then identify Sym(n) with a quasi-projective subset of the complex projective variety 
P^. 

Definition 3.7. The algebraic closure of the set 

{rowsp [F In] \ F E Sym(n)} 

is called the complex Lagrangian Grassmann manifold. It will be denoted by LG(n). 

It is well known that lLG{n) is a smooth projective variety of of dimension ("^^), the 
dimension of Sym(r;,). Note that every element in hG{n) can be simply represented by a 
subspace of the form rowsp [Fi F2], where ^1(^2)* is a symmetric matrix, i.e. Fi{F2Y = 
-^2(-^i)*- The elements of hG{n) are thus exactly the Lagrangian subspaces of C^*^. The 
subspace rowsp [Fi F2] coincides with the subspace rowsp [5* /„] associated with an element 
5" of Sym(?T,) if and only if F2 is invertible. Moreover, then 5* = {F2)~^Fi. When F2 is 
singular one can still define a characteristic polynomial through 



(p{s) := det 



D{s) N{s) 
Fi F2 



(3.5) 



Note that in the Hamiltonian case, (f{s) is necessarily even, i.e. then 

<^{s) = v?(-s). 

Let fi, i = 0, . . . , N he the Pliicker coordinates of rowsp [Fi F2]. In terms of the Pliicker 
coordinates the characteristic equation can then be written as: 



det 



D{s) N{s' ^ 
Fi F, 



2 



$^P.(5)/., (3.6) 



j=0 



where Pi{s) is the cofactor of fi in the determinant (j3.6p . 
Let Z C P^ be the linear subspace defined by 



N 



Z = {^gP^|5^j9,(s)z, = 0}. (3.7) 



j=0 



Following [ini I2SI ISni we identify a closed loop characteristic polynomial ip{s) with a 
point in P"^. In analogy to the situation of the static pole placement problem considered 
in (compare also with Section 5]) one has a well defined characteristic map 



X : LG'(n) ~Z — > P 



s 



rowsp [Fi F2] I — > E*=o fiP 



{3.t 
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in the complex symmetric case and 



X' : hG{n) - Z 
rowsp [Fi F2] 



p<5/2 

even part of Xlio fiPii^)- 



(3.9) 



in the Hamiltonian case. In the latter case the reduction in dimension of the projective space 
arises due to the evenness of the closed loop characteristic polynomial, so that in the second 
map only the coefficients of the even terms of fiPii^) ^'^ appear. 

Recall the notion of degree of a variety fHl Chapter I, §7] and the notion of a central 
projection (see ^W*, Chapter I, §4]). The geometric properties of the map x ^^e as follows: 

Theorem 3.8. The maps Xi x' define central projections. In particular if Z {~\ LG'(n) = 
and dimLG(?T,) = ("2^) = 5 then x is surjective, and there are degLG(n) many pre- 
images (counted with multiplicity) for each point in P*^, where degLG'(ra) is the degree of 
the Lagrangian manifold l^G{n) in P^. Similarly, z/dimLG(n) = ("'^^) = (5/2 then x' is 
surjective with exactly degLG(n) many pre-image points in each fiber. 

Proof. By definition (see e.g. ISSIEOI) X represents a central projection of 'LG{n) from the 
center Z to P*^. When Z fl LG(n) = and dimLG'(n) = ("^^) = 5 then x is a finite 
morphism [301 Chapter I, §5, Theorem 7] and onto of degree degLG(n) |23 Corollary 5.6] 
Similarly for x'. □ 

The set Z fl LG'(n) is sometimes referred to as the base locus. The interesting part of 
the theorem occurs when the base locus Z fl LG(n) = since in this situation very specific 
information on the number of solutions is provided. \i Z H LG(n) = and ("^^) = ^ (or 
n{n + 1) = 5) then one says that x (or x') describes a finite morphism from the projective 
variety LG(n) onto the projective space P'' (or P'^/^). 

This last situation is most desirable and this motivates the following definition. 

Definition 3.9. A particular symmetric transfer function G{s) is called nondegenerate if 
Z n LG(?2) = 0. A system which is not nondegenerate will be called degenerate. 

In terms of matrices a symmetric transfer function G{s) = D{s)~^N{s) is degenerate as 
soon as there is a Lagrangian subspace rowsp [Fi F2] G LG(n), such that 



In a slightly more geometric language this means that the Hermann- Mart in curve 
defined by rowsp [D{s) N{s)] is fully contained in a Lagrangian hyper-plane defined by 
rowsp [Fl F2] . In the study of the static pole placement problem [3] and the dynamic pole 
placement problem ^] definitions analogous to Definition 13.91 plaved an important role. 

The next lemmas give specific information under what conditions Z fl hG{n) = 0, i.e. 
under what conditions a symmetric transfer function is nondegenerate. Similar results were 
crucial in proving the pole placement results in IT^ I77j. 

Lemma 3.10. If 6 < ("^^) = dimLG(n) then every nx n symmetric transfer function 
of McMillan degree 6 is degenerate. Similarly, any n x n Hamiltonian transfer function of 
McMillan degree S is degenerate, if S < n{n + 1) = 2dimLG(n) 



det 



N{s) 
F2 



0. 
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Proof. dimZ >N — 6 — lasZis defined by (5 + 1 linear equations (5/2 + 1 many in 
the Hamiltonian case). If dimLG(n) > 6 (or dimLG'(?2) > 6/2 in tlie Hamiltonian case), 
tfien Z nLG(n) is nonempty by tlie (projective) dimension tlieorem (see e.g. jT8| Cliapter I, 
Tlieorem 7.2]). □ 

Lemma 3.11. If 5 = ("^^) = dimLG'(n) (or 6 = n{n + 1)), then a generic set of n x n 
symmetric (or Hamiltonian) transfer function of McMillan degree 6 is nondegenerate. 

Proof. Let Q be tlie set of all n x n symmetric transfer functions of McMillan degree n. Q 
can be given the structure of a quasi-projective variety. For this recall the definition of the 
projective variety Kf^ ^ introduced in j27j and which compactifies the set of all n x n transfer 
functions of McMillan degree 6. An element (Hermann- Mart in curve) rowsp [D{s) N{s)] G 
K^n describes an element of Q as soon as degdetD(s) = 6 and D{s)N{sy = N{s)D{sY. 
The last condition translates into some linear conditions to be satisfied among the Pliicker 
coordinates of Kf^n- The resulting sub-variety of K^n constitutes a natural compactification 
of Q and Q itself is a quasi-projective subset. 
Consider now the coincidence set 



S := \ {D{s)-^N{s); ^1,^2) G Q x LG(n) | det 



D{s) 
Fi 



N{s) 
F2 







Since hG{n) is projective the projection onto Q is an algebraic set by the main theorem of 
elimination theory (see e.g. ^3)- The set of nondegenerate systems forms therefore a Zariski 
open subset of Q. We have shown the result if we can exhibit one n x n transfer function of 
McMillan degree ("^^) which is nondegenerate. The next lemma gives such an example and 
the claim therefore follows. Note that the previous arguments run completely similar for the 
Hamiltonian case and it therefore remains to construct one example as well. However, the 
symmetric Hamiltonian transfer function G'(s^) does exactly the job. 



Lemma 3.12. The symmetric transfer function G{s) 



ate. 



□ 



is nondegener- 



Proof. First it is clear that G{s) has McMillan degree 6 = {^2^) and that 



[Dis) Nis)] 



forms a left coprime factorization of G{s). Let 



R :-- 
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and assume by contradiction that G{s) is degenerate. It therefore exists rowsp [Fi F2] G 
L,G{n), such that 



= det 



D{s) 
Fi 



N{s) 
F2 



det 



D{s) 
Fi 



N{s)R 
F2R 



(3.10) 



Let 5* G Gin be the matrix which transforms the n x 2n matrix [Fi F2R] into row reduced 
echelon form, i.e. 



{SF2R)] 



Let 



* •■■ * 1 ■■• ■■• ••■ ■•■ 

* ■■■ * * ■■■ * 1 ■■■ ■■■ ■•■ 



* * 



* 



* 1 



=: [Fi F2 



(3.11) 



ii < . . . ik < n < ik+i < ... in < "in 

be the pivot indices. We claim that the first k pivot indices determine the last n — k pivot 
indices uniquely. For this let < . . . < i„„fc be the complementary indices of the indices 
{«!,... ik} inside the set {1, . . . , n}. Then we claim that: 



^k+l 



2n - in-k + 1 



in = 2n — ii + 1. 

Indeed, if this is not the case then it follows that FiR{F2y cannot be symmetric for any 
choice of values in the row reduced echelon form ()3.1H1 . On the other hand the matrix 
FiR{F2y has to be symmetric since by assumption ^1(^2)* is symmetric. 

The indices zi, . . . , z„ describe the maximal Pliicker coordinate (with regard of the Bruhat 
order) of rowsp [Fi -F2-R] which is nonzero and the correponding cofactor of [D{s) N{s)R\ 
is computed as ±s°, where a = ^"jf In general there are other fullsize minors (Pliicker 
coordintes) of [D{s) N{s)R\ which have the form ±s". All other Pliicker coordinates with 
this value are however not comparable with regard to the Bruhat order and since ii, . . . ,in was 
the maximal nonzero Pliicker coordinate of rowsp [Fi F2R\ it follows that the determinant 
expansion in ()3.10|) cannot be zero. This is a contradiction and it follows that G{s) is 
nondegenerate. □ 

Remark 3.13. For the static pole placement problem Brockett and Byrnes ^ showed that 
the osculating normal curve 



rowsp 



1 
... 



2s ^n.+p-l^^ 



j-m+p— 1 
m+p-l\ rn+p-2 







(m+p~l\ ^p 
\ m-1 J 



G Grass(m, 



'"m+p\ 



is nondegenerate. Also in this situation the Pliicker coordinates have the simple form ±s^, 
where /3 = YlT=i ~ ^ there are no two Pliicker coordinates which are comparable in 
the Bruhat order and give rise to the same monomial s^. 
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We have now all pieces together in order to prove the main result. 

Proof of Theorem \3.(A W.l.o.g. we focus on the case of symmetric transfer functions. The ar- 
guments for the Hamiltonian case run completely similar. Note however, that the closed loop 
characteristic polynomial of a Hamiltonian system is always an even polynomial. Therefore 
our definition of generic pole-assignability for Hamiltonian systems restricts to the space of 
even polynomials. Since the dimension of the space of even monic polynomials of degree 5 is 
5/2, the appropriate condition for Hamiltonian systems is 5/2 < ("^^). With the comments 
in mind we return to the proof for symmetric transfer functions. 

When 5 > ("^^) then a simple dimension argument shows that the image of the charac- 
teristic map X described in fl3.8p has dimension at most ("^^) and therefore there is a Zariski 
open set in P'' not in the image of x- 

When 6 = ("^^) then Lemmas 13.111 and 13.121 show that there is a generic set of n x n 
symmetric transfer functions of McMillan degree 6 which are nondegenerate. The character- 
istic map ()3.8|1 has therefore no base locus and every point in the image of x has deg]LG{n) 
pre-image points when counted with multiplicities. The degree of the variety hG{n) was 
recently computed by Totaro and it resulted in the number ()3.3|) . 

A priori the geometric formulation only predicts degLG(n) many solutions inside lLG{n) 
and it is not clear if all these solutions correspond to regular feedback laws of the form 
u = —Fy + V. If G{s) is a strictly proper symmetric transfer function then this is indeed the 
case and the same argument applies as in [3]. □ 
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